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Abstract 

This paper proposes a new approach to torce 
balance structural optimization featuring a 
computational experiment design. Currently, this 
multi'dimensional design process requires the designer 
to perform a simplification by executing parameter 
studies on a small subset of design variables. This one- 
factor-at-a-time approach varies a single variable while 
holding all others at a constant level. Consequently, 
subtle interactions among the design variables, which 
can be exploited to achieve the design objectives, are 
undetected. The proposed method combines Modern 
Design of Experiments techniques to direct the 
exploration of the multi-dimensional design space, and 
a finite element analysis code to generate the 
experimental data. To efficiently search for an 
optimum combination ot design variables and minimize 
the computational resources, a sequential design 
strategy was employed. Experimental results Irom the 
optimization ot a non-traditional force balance 
measurement section are presented. An approach to 
overcome the unique problems associated with the 
simultaneous optimization of multiple response criteria 
is described. A quantitative single-point design 
procedure that reflects the designer’s subjective 
impression of the relative importance of various design 
objectives, and a graphical multi-response optimization 
procedure that provides further insights into available 
tradeoffs among competing design objectives are 
illustrated. The proposed method enhances the intuition 
and experience of the designer by providing new 
perspectives on the relationships between the design 
variables and the competing design objectives 
providing a systematic foundation for advancements in 
structural design. 
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Introduction 

Direct force and moment measurement ot 
aerodynamic loads is lundamental to wind tunnel 
testing. Typically these measurements are made with 
an instrument known as a force balance, which is 
internally mounted in a scaled wind tunnel model (see 
Figure 1 ). The force balance is a complex structural 
spring element, which provides high-precision 
measurements ot the aerodynamic loads exerted on the 
wind tunnel model, by measuring strain within its 
flexural elements. This measured strain is then 
converted into an electrical signal. Electrically 
measured strain as a tunction ol an externally applied 
load forms the basic concept of force balance 
measurements. 



Figure 1. Internally mounted force balance. 


The balance flexural elements are designed such 
that the magnitude ot the strain, at strategic locations, is 
the same for the individual application ot each 
component of load, even though the magnitudes of the 
components are not equal. Therefore, the flexural 
element stiffness must be simultaneously tuned with 
respect to the component loads in six degrees ot 
freedom. An adequate balance design is achieved 
through a combination of the structural design ot 
multiple flexural elements, and the strategic location of 
the strain gages. Auxiliary factors must also be 
considered by the balance design engineer including the 
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complexity of fabrication, which is directly 
proportional to both the production time and cost of the 
force transducer. Ultimately, the physical design is 
constrained by the available space within the wind 
tunnel model. 

Currently, structural optimization of the force 
balance relies heavily on the experience and intuition of 
the balance designer. Software based analysis and 
simulation tools have been introduced in recent years to 
aid in this process', but the decisions regarding design 
variables are solely driven by the balance designer. 
These design variables (independent variables) include 
the measuring section dimensions, shape, location, and 
strain gage placement. 

The quality of the design is evaluated based on 
multiple competing response metrics (dependent 
variables) such as: achieving an acceptable sensitivity 
of the measurements, maximizing the stiffness, 
minimizing the combined stress intensity, and 
minimizing the fabrication complexity. There is 
typically on the order of 20 independent variables, and 
10 dependent variables. The goal of this multi- 
dimensional design process is to find an optimum 
combination of the design variables that creates the best 
overall performance of the response metrics within the 
constraints. 

It is very difficult to reconcile the effects of such a 
large number ot design variables on so many comp>eling 
response metrics by conventional one-factor-at-a-time 
(OFAT) methods. In practice the balance engineer 
typically uses his intuition to effect a simplification of 
the problem. This simplification is often achieved by 
considering only a relatively small subset of the design 
variables, and focusing primarily on a few of the 
response metrics believed to be most critical. Even in 
this case, OFAT methods make it difficult to resolve 
potentially important interactions among the design 
variables. 

Past research in this area has identified the 
structural optimization of a force balance as an 
emerging area of force balance development. Zhai" and 
others identified the need for optimizing the structure of 
the balance to achieve the desired high sensitivity, low 
interactions, and high stiffness. An empirical 
algorithm, which simulates biological growth, was 
combined with classical finite element analysis to 
develop new llexure shapes. Zhai^ and Parker'^-^ 
highlight the benefits of a stiffer balance, namely that it 
will raise the natural frequency of the model, balance, 
sting system thereby reducing undesirable model 
vibration. The desirability of high sensitivity, or 
flexibility, and high stiffness emphasizes the competing 
criteria involved in the structural optimization. 
Parametric studies of beam dimensions were performed 


using an OFAT approach to develop optimum flexural 
geometry.^ An important finding of this effort was that 
only minor improvements are possible without 
investigating new structures and shapes. 

Iterative gradient-based algorithms were applied in 
Reference 6. In this effort, Hou and Twu identified the 
ability to automate the design of balance flexural 
sections. The main goal of this work was to maximize 
the sensitivity, while minimizing the stress intensity. 
Once again, the competing nature of the criteria is 
highlighted by the desire to have high stress levels in 
the local vicinity of the strain gage, while minimizing 
the stress globally over the flexural element. 

Approach 

This paper proposes the use of a computational 
experimental design using Modern Design of 
Experiments (MDOE) techniques combined with a 
state-of-the-art finite-element analysis code to generate 
the experimental data. The goal of the proposed 
method is to efficiently move through the multi- 
dimensional design space incorporating the competing 
criteria for the responses. 

The proposed technique augments the intuition of 
the designer by providing graphical insight to the 
relationship between the design variables and the 
response metrics. This insight differs from other 
structural optimization methods that could be used such 
as: genetic algorithms, neural networks, and gradient 
search techniques. These methods search the design 
space to find a maximum response location, but many 
times do not provide vital information on why a 
particular path through the design space was chosen. In 
contrast, the current effort builds mathematical models 
of a portion of the design space providing an 
representation of its landscape. These mtxiels provide 
the balance designer with the relative influence of the 
design variables and valuable information about the 
subtle interactions among them. Other optimization 
methods commonly provide an answer without an 
explanation, while the use of experimental design 
techniques offers both an answer, an explanation, and 
new insights to the balance design engineer. 

A sequential design strategy was employed in order 
to minimize the required number of computational runs, 
thereby increasing computational efficiency. This 
sequential strategy relies on lower order mathematical 
mcxlels initially, that account for the main effects of the 
design variables and the interactions among them. 
These simple, initial models serve as a relatively low- 
cost means of identifying efficient directions for further 
experimentation. This strategy then involves an 
increase in the order of fitted models in the more 
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interesting regions of response extrema (maximum or 
minimum responses.) Once an adequate model has 
been generated, an optimization is performed to 
determine the combination of design variables that 
maximizes a weighted desirability function that 
combines multiple response metrics. 

Experimental Results 

This section details an example of the optimization 
of a force balance measurement section using the 
proposed technique. A delinition of the design and 
response variables is provided followed by a discussion 
of the experimental design and analysis. The final 
process of using the derived mathematical model to 
perform an optimization is presented. 


Geometry Description 


A typical six-component force balance used at 
NASA Langley Research Center consists of three 
measurement sections as shown in figure 2. These 
measurement sections are the areas of concentration for 
the structural design. The center section is designed to 
measure axial force, or drag. The two sections on either 
side of the axial section are designed to measure the 
other five components of load and are referred to as the 
cage sections. The term cage section comes from the 
classical geometry of multiple rectangular flexures used 
in this section creating a cage type structure. A parallel 
system of rectangular flexures is used because it is 
relatively easy to analyze, but does not necessarily 
represent optimum geometry. In fact the sharp corners 
of the rectangular flexures create locations of stress 
concentration, especially due to torsional loads. On 
either side of the cage section are solid elements of the 
balance structure, which are referred to as bulkheads. 
These bulkheads are considered rigid compared to the 
flexural elements of the measurement sections. The 
balance moment center, which is a reference location 


for the torque components, is typically defined to be at 
the center of the axial measurement section, and 
therefore the cage sections are symmetric about the 



Figure 2. Typical force balance measurement sections. 


An advanced geometry cage section was selected 
ibr optimization. As previously mentioned, minimizing 
the number of flexures that comprise the cage section is 
desirable from a manufacturing point of view; therefore 
a single flexural element was considered. It is believed 
that a hollow elliptical cross section would be optimum. 
For this experiment, the elliptical shape has been 
approximated by an octagon with a hole in the center. 
This approximation is due to manufacturing difficulty 
of an elliptical beam. The design variables control the 
volume, shape, contour, and location of the cage 
section. 

A generic octagonal cross section is shown in 
figure 3a. The four linear dimensions of the cross 
section, the diameter of the center hole, and the two 
radii that provide a smooth transition between the sides 
of the octagon represent seven of the design variables. 
A side view of the flexural element is provided in figure 
3b. In this figure, the location of the cage section with 
respect to the balance moment center, the length of the 
flexure, and the symmetrical radii at the ends of the 
flexure represent three additional design variables. A 
summary of these ten design variables is provided in 
Table I. 



(a) Octagonal cross section. 
Figure 3. Design variables. 



(b) Side view of the flexure. 
Figure 3. concluded. 
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Table 1. Summary of the design variables 


Design 

Variable 

Description 

A 

location of cage section 

B 

minor width 

C 

minor height 

D 

overall width 

E 

overall heisht 

F 

length of flexure 

G 

hole diameter 

H 

radii at ends of flexure 

J 

longitudinal radii along lop of flexure 

K 

longitudinal radii along side of flexure I 


The response variables determine the quality of the 
structural design of the force balance. As previously 
mentioned, the goal is to maximize the sensitivity while 
minimizing the structural deflection and overall stress 
intensity. The four measurement sensitivities chosen 
for this experiment were normal force, pitching 
moment, yawing moment, and side force as shown in 
figure 1. The strain gages used to measure these 
components are located at the longitudinal centerline of 
the flexure on the surfaces with the maximum distance 
Irom the neutral axis with respect to the component of 
load. The strain gage locations are shown in figure 4. 
The balance utilizes an electrical network of four strain 
gages configured in a Wheatstone bridge arrangement. 
In this experiment a single location was used to 
determine the sensitivity due to the symmetric 
placement of the four gages in the bridge. Furthermore, 
this component measurement strategy uses the same 
gage to measure normal force as well as pitching 
moment. Normal force is measured by adding the 
signals from the two symmetric cage sections (forward 
and aft of the balance moment center) and pitching 
moment is measured by subtracting the signals. The 
yawing moment and side force are measured in a 
similar arrangement. 



Figure 4. Response measurement locations. 


The rolling moment component, which would 
typically be measured in the cage section, was not 
included in this initial investigation for the purpose of 
simplification. The location of the strain gage for 
rolling moment could be anywhere on the flexure 
depending on the cross section. This would increase 
the number of design variables to at least thirteen, and 
add considerable complexity to the automated aspects 
ol the finite element solution and results interpretation. 
Including rolling moment in future work would provide 
a significant advantage over current methods. 

Additional response variables were included that 
were used to minimize the structural deflection and the 
stress intensity. The structural deflection is measured at 
a location forward of the cage section as shown in 
figure 4. There are two response variables for the 
deflection, one in the normal force plane and one in the 
side force plane. It is desirable to minimize both of 
these simultaneously. The stress intensity was 
calculated by superposition of the maximum von Mises 
stress due to the application of the four components of 
load under study (normal, pitch, yaw, and side). The 
location of the maximum stress was not considered for 
each flexure configuration, but intuitively it will be 
located at the transition between the cage section and 
the bulkheads. A summary of the respon.se variables is 
provided in Table II. 


Table II. Summary of the response variables 


Response 

Variable 

Description 

NF_out 

normal force sensitivity 

PM_out 

pitching moment sensitivity 

YM_oul 

yawing moment sensitivity 

SF^out 

side force sensitivity 

zdisp 

displacement in the normal plane 

ydisp 

displacement in the side plane 

max_stress 

maximum Von-mises stress 


To summarize the experiment parameters, there are 
ten design variables (independent variables) and seven 
response variables (dependent variables). We want to 
obtain a mathematical model that enables the prediction 
of the values for the response variables from a given 
combination of the design variables. With this model 
we are able to search through the design space to find 
an optimum combination of design variables to achieve 
the combined response objectives. The experimental 
data required to generate the model was obtained using 
a finite element analysis package described in the next 
section. 
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Solid Geometry Modeling and 
Finite Element Analysis Description 


changes in the flexure geometry without manual 
intervention. 


A parametric solid model was used to generate the 
geometry for the finite element analyses. The 
representation of the force balance was simplified to 
minimize computational resources. The force balance 
was longitudinally cut in half at the balance moment 
center. All of the features except those necessary to 
model the flexure boundary conditions, apply the loads, 
and properly constrain the system were removed. The 
simplified geometry is shown in figure 5. A table of 
parameters that were specitied by the experimental 
design controlled the geometry. Each combination of 
design variables was verified to ensure that it did not 
violate the geometrical constraints. This process of 
verification will be described in more detail in the next 
section describing the experimental design. 


Constrained surface - 




-Cylindrical surface used 
to apply the loads 


Figure 5. Parametric solid geometry. 

The finite element analysis (FEA) package used for 
this experiment is based on p-element technology. 
elements rely on a relatively coarse discretization ot the 
geometry using solid tetrahedral elements. Differing 
from classical finite element analysis codes, the transfer 
functions through these elements are high-order 
polynomial functions. This type of FEA has the 
desirable property of self-refining its ability to represent 
steep stress intensity gradients found in geometry 
transitions. In a classical FEA system the analyst must 
manually refine the mesh in these areas of stress 
concentrations. In a p-element code, the elements are 
unchanged but the order of the polynomial transfer 
function is iteratively increased up to order nine. This 
is done internally during the iterative solution process 
and does not require manual intervention. This type of 
system was well suited for the task ot generating the 
experimental data since it was robust in handling large 


Experiment Design and Analysis 

This section describes the construction and analysis 
of the experimental designs used to generate the 
mathematical models of the design space. Factorial 
experimental designs were used because they provide 
an efficient means to determine the influence and 
interaction among the design variables on the 
responses. Even though this is a computational 
experiment, the data is not without cost. Computational 
resources and manual manipulation of the results are 
required to obtain the data, and therefore the data 
volume should be minimized- Factorial experiment 
designs are well suited for this application. 

In a classical one-factor-at-a-time approach, a 
single design variable would be varied while holding all 
other variables at a constant level. While this approach 
is intuitive to the designer, it is inefficient and more 
importantly could miss significant interactions among 
the design variables. With the OFAT approach, it is 
assumed that the effect of a single variable on the 
responses is the same regardless of the level ot the other 
design variables. In real world situations, this is rarely 
the case. Factorial designs vary all of the variables 
simultaneously in a specified manner and enable the 
determination of the same main effects derived from an 
OFAT approach as well as significant interactions 
between the variables. 

A two-level factorial design was used initially. 
These designs are useful in a broad array of 
applications but they are limited in one important 
respect: The fact that they feature only two levels of 
each design variable means that response models 
developed from data acquired with such designs can 
only accommodate first-order and mixed first-order 
terms. For example, in a two-variable, two-level 
factorial design with factors A and B, it is possible to 
develop response models featuring the linear A and B 
terms and a second-order interaction term, AB. The 
pure quadratic terms (A" and B') cannot be quantified, 
because this would require at least three levels of the 
design variables. 

A useful extension to the two-level factorial 
design, accommodating full second-order response 
models, is the Box-Wilson, or Central Composite 
Design (CCD). A CCD for two variables is illustrated 
in figure 6. The four data points comprising the square 
are the two-level factorial design points and are referred 
to as “corner points'’. They define all combinations ot 
the high and low levels of the two variables. There are 
also a number of “center points' in this design that are 
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replicates of the variable levels at the design center. 
Finally, there are four “axial points", also referred to as 
“star points", arranged on the axes of the variable 
coordinate system. A CCD features five levels of each 
design variable. 



Factor A 

Figure 6. Central Composite Design in two variables. 

Although the description of the CCD has been 
given in a two-variable case, the design extends readily 
to any number of variables. While the details of the 
CCD are beyond the scope of this paper, many 
advantages flow from this design space geometry. The 
reader is referred to References 7, 8, and 9 for more 
details. 

The first step in constructing the experimental 
design was to select the range of the design variables. 
The selection of the range proved to be a challenging 
aspect of this experiment. This was due to the need of 
choosing a wide enough range while simultaneously 
ensuring that the geometrical constraints were fulfilled. 
In general, to produce a useful model from an 
experimental design, the design variables must be 
varied over a large enough range to effect a meaningful 
change in the responses. For the initial design, the 
range was chosen to be as broad as po.ssible. In 
subsequent designs, as we moved toward a response 
extrema, the range was narrowed. 

Furthermore, the combination of design variables 
was constrained due to the physical geometry that they 
represent. An arbitrary combination of variable levels 
could violate the physical geometry. For example, the 
hole diameter must always be less than the overall 
width or height of the flexure. If the diameter were 
larger than one of these dimensions, it would cut 
through the sides of the flexure and thereby produce 
unacceptable geometry. Prior to executing the design, 
all of the combinations of the design variables were 
verilied to ensure that the geometry was acceptable. If 
a particular combination was not acceptable then it was 
slightly modified instead of adjusting the range and 
generating a new design. Tailoring the design to meet 


the geometrical constraints proved to be tedious and 
somewhat limiting during the experimentation. The 
ability to mathematically define the physical constraints 
and build them into the experimental design is an area 
of future research. 

Three experimental designs were performed to 
determine a final flexure design. The first experiment 
involved the variables A, D, E, and F. These variables 
control the volume and location of the flexure. The 
goal of this first experiment was to define the center of 
the design in the search space. The design center is 
represented by a combination of variables that is near a 
combined response extrema. For this first design, the 
flexure was rectangular in shape. A four-variable full 
factorial design consisting of 16 runs was performed. 
The 16 geometries represented by this design are shown 
in figure 7 (the forward portion of the balance has been 
removed for clarity). Mathematical models consisting 
of the main effects and two-way interactions were 
derived. Using these models, the design center was 
approximated for these four variables. 

The second design included the other six design 
variables that control the shap>e and contour of the 
flexure. Depending on the level of the design variables, 
the shape of the flexure was able to transition from a 
rectangle to a thin-walled octagonal cross section. A 
highly fractionated factorial design was executed to 
minimize computational resources. With 10 design 
variables, a full factorial would require 1024 runs. Not 
only is this an excessive number of runs, but also it is 
not necessary to derive an adequate model. A full 
factorial design in 10 variables allows for the resolution 
of all possible interactions including the main effects, 
two-way, three-way, and four-way, up to ten-way 
interactions. While it is possible that there are 
significant two-way and three-way interactions, it is 
unlikely that ten-way interactions would make a 
significant contribution to the model. 

Therefore, a fractionated design enabled fewer 
experimental runs while giving up the capability to 
determine some multiple factor interactions. The 
consequence of a fractionated design is that some lower 
order terms will be aliased with higher order terms, 
thereby making it difficult to unambiguously resolve 
the lower order terms. A 1/64**" fraction of the full 
factorial was executed requiring only 16 combinations 
of the 10 design variables. The geometries represented 
by these 16 combinations are shown in figure 8. The 
mathematical models derived from these design were 
then used to determine combinations of the design 
variables that would fulfill the response criteria. The 
analysis of this design revealed that it contained a 
combined response extrema. 
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Therefore, in the third and final experiment, the 
fractionated factorial design in 10 variables was 
augmented to include axial and center points. This 
augmentation created a hybrid form of a central 
composite design. This design combined the 16 runs 
from the 1/64“’ fraction with 20 axial points and 10 
center points resulting in 46 total points. There were 
now five levels of each of the design variables. A 
single center point geometry was used and then 


simulated random errors were added to generate ten 
replicates. The inclusion of axial and center points 
enabled the determination of curvature in the math 
model. Therefore, the use of higher order terms 
increased the precision of the prediction and was 
important for the optimization of the variable levels. 
The geometries of the 21 additional points are shown in 
figure 9. 



Desirability Function 

One complication in an optimization problem that 
features multiple response criteria is that the 
combination of independent variables that optimizes the 
design with respect to one response is not generally the 
combination that optimizes it for another. Derringer 
and Suich'“ developed the concept of a “desirability 
function"' to cope with such situations. The desirability 
function is a single quantitative response index ranging 
from 0 to 1 that approaches I more closely as various 
criteria specified by the designer are more closely met. 
Desirability functions are computed for each response 
individually. Calculating the geometric mean of all the 
individual functions then produces a single index. The 
desirability function can identify a solution that consists 


of some combination of independent variables that may 
not perfectly optimize the performance of the design 
with respect to any of the individual responses of 
interest, but which provides desirable performance for 
all of the responses. 

For example, assume that it is unacceptable if a 
response level falls below some value. A, and that the 
largest possible value between A and some goal value, 
B, is desired. Assume further that any response greater 
than or equal to B is equally acceptable. That is, we 
assume that the designer feels that a given response has 
to exceed the A level to be minimally acceptable and 
the larger the better after that, with all expectations met 
il the response level is B or more. The desirability 
function for this response variable is created by a 
transformation mapping the range from A to B into a 
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range from 0 to 1. Response values less than A are 
assigned a desirability value of 0, and those greater than 
B are assigned a value of 1. If the desire is to achieve 
the largest possible value of the response, the value of 
“B" can be set to the largest value of the response that 
can be achieved within the full range of independent 
variable settings. 

Similarly, if response levels greater than some 
value, B, are unacceptable, and the smallest possible 
value between B and some minimal goal value. A, is 
desired, and if any response as small or smaller than A 
is equally acceptable, the desirability function is created 
by a transformation mapping the range from A to B into 
a range from 1 to 0. Response values greater than B are 
assigned a desirability value of 0, and those less than A 
are assigned a value ot 1 . It the desire is to achieve the 
smallest possible value of the response, the value of 
“A ’ can be set to the smallest value ot the response that 
can be achieved within the full range of independent 
variable settings. 

It may simply be desirable tor a particular response 
to lie with some range of values between A and B. In 
that case, the desirability (unction is given a value of 1 
for responses between A and B and 0 tor responses 
outside this range. 

It is often desirable that a particular response lies as 
close to a specific target value as possible. In this case, 
the desirability tunction is designed to have a value of 1 
at the target value and to approach zero as the response 
varies from the target in either direction, becoming zero 
when the response is unacceptably distant from the 
target. 

It is not necessary for the desirability function of a 
specific response to vary linearly from 0 to 1. A 
weighting constant, w, can be applied that imparts to 
the desirability function the designer's subjective 
evaluation of its role in the overall optimization. For 
example, the desirability tunction for a response to be 
maximized in the range ot A to B can be constructed to 
have a value of 0 it the response is less than A, 1 it the 
response is greater than B, and tor responses between A 
and B, the desirability function can be described as 
follows: 


B-A 


where d, is the desirability function for the A' response, 
and vr is a weighting constant that is smaller when 
maximizing the response is relatively less important 
and larger when maximizing the response is relatively 
more important. Figure 10 shows how this function 
changes as the response varies between A and B. for 


various values of vv. Similar weightings can be applied 
for other objectives besides maximizing the response. 



Figure 10. Desirability function 
for a response to be maximized. 


After desirability values are computed for each 
response variable, they arc combined into a single 
desirability index, D, by calculating their geometric 
mean. Further refinements to the weighting can be 
applied by assigning a range ol numbers (from I to 5, 
say) to the importance of optimizing each response 
variable. The final desirability index then is computed 
as follows: 


0= IK 





where v, is a number indicating the relative importance 
of the response, which might typically be an integer 
in the range ot I to 5, with 5 indicating the greatest 
importance and 1 indicating the least. 

In the present study, desirability functions were 
developed for the criteria that sensitivities tor normal 
(bree, side force, yaw moment, and pitch moment all 
exceed 1000 microvolts per volt and be as large as 
possible, while the deflection ot the balance in the z and 
y axis be minimized, and in any case less than 0.002 
inches in the z axis and 0.(X)5 inches in the y axis. The 
maximum stress intensity was specitied to be less than 
75,000 pounds per square inch, and as small as 
possible. Equal importance was assigned to all 
response criteria (v^l for all /) and linear weightings 
were applied (h-I.O for all responses). Combinations 
of the independent variables that maximized the grand 
desirability function, D, were chosen as the starting 
point for further refinements applied using graphical 
optimization, as outlined in the next section. 
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Graphical Optimization 

The desirability functions are based on the 
designer’s subjective impression of the relative 
importance of various design objectives, and can best 
be thought of as a tool to achieve an approximately 
optimum single-point design. Further refinements on 
this initial design are possible using graphical 
optimization techniques that offer the designer 
additional insights into the tradeoffs among competing 
optimization criteria. These methods result in a range 
of independent variables throughout which all 
responses satisfy certain minimal criteria. 

in this study, an initial flexure design was 
determined using desirability analysis as described 
above. Special computer software was then used to 
perform graphical optimization by exploiting the fact 
that extrema in the various response variables have 
associated with them closed contours of constant 
response value. Even if such an extremum lies beyond 
the range of independent variables considered, there 
will often be open contours that enclose some region of 
the design space bounded on one or more sides by a 
design space boundary. 

It is possible to define contours for each response 
variable that correspond to some particular optimization 
criterion. For example, the contours can be defined that 
enclose all combinations of the independent variables 
lor which a given response exceeds some minimum 
acceptable value or is less than some maximum 
acceptable value. Pairs of contours can be selected to 
bound combinations of independent variables for which 
the response is within some prescribed range. 
Overlaying such contours for different response criteria 
reveal at a glance whether there is a non-null 
intersection of all the criterion regions. If not, then it is 
possible to see which criteria require the least relaxation 
to achieve such an overlap. It is also possible to see 
that a relatively inconsequential relaxation of one or 
more criteria might dramatically increase the number of 
combinations of independent variables that satisfy all 
constraints. It could be quite beneficial to have a fairly 
wide range of independent variables that meet all 
response criteria, since certain combinations may be 
easier to implement than others or less expensive. 
Figure II is representative of a graphical optimization. 
The area highlighted in gray indicates combinations of 
the independent variables that satisfy all constraints. 
This region is known colloquially as a "sweet spot”, 
and graphical optimization is therefore often described 
informally as "sweet spot analysis”. 

Once a region is defined by graphical optimization 
within which all response criteria are satisfied, a 
specific design point can be selected within this region 


to maximize or minimize some other factor of interest. 
For example, not all balance geometries that satisfy 
performance goals will be equally easy to manufacture. 
In such a case, the combination of geometric 
parameters that minimizes the manufacturing 
difficulties can be selected from within the region that 
satisfies all performance requirements. Developing 
such manufacturing cost functions is a next logical 
extension of the optimization procedures described in 
this paper. 



Figure 11. Graphical optimization. 

In this case figure 1 1 illustrates design variable F 
(length of the flexure) plotted versus variable A 
(location of the flexure). To overcome the limitation of 
not being able to plot all of the design variables 
simultaneously, which would require a ten dimensional 
plot, two are selected and all other variables are set to 
the constant levels. The contour lines represent the 
response criteria previously specified. The response 
boundaries that border the sweet spot are the normal 
force sensitivity, side force sensitivity, and the 
maximum stress intensity. All other response criteria 
are some distance from the boundaries of the sweet spot 
indicating that they are not limiting metrics for this 
particular combination of design variables. Any 
combination of variable A and F lying within the sweet 
spot will satisfy all of the specified criteria. 

Based on this graphical optimization, the final 
flexure geometry was selected and is illustrated in 
figure 12 (only the flexure is shown for clarity). A 
combination of design variables in the upper left corner 
of the sweet spot was selected because that area 
represents the minimum stress intensity, while 
satisfying the sensitivity requirements. A finite element 
analysis was performed on the resulting geometry and 
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was compared to the prediction oi the mathematical 
models derived from the experimental design. The 
comparison is provided in Table III. All ot the 
predicted values fall within the 99% prediction inicrvals 
of the model, except for the prediction of the 
displacement in the y-axis, which is slightly outside the 
interval. This level ol agreement between the 
prediction and the actual FEA solution demonstrates the 
usefulness ol the mathematical models lor searching 
through the design space. 



Figure 12. Final llexure geometry. 


Table III. Predicted responses and FEA results. 


Response 

Variables 

Response Model 
Prediction 

FEA solution 

NF out 

1004 

982 uVA' 

PM out 

1080 

1052 uVA-' 

YM_out 

1091 |iVA/ 

1152nVA^ 

SF_out 

1008 uVA^ 1 

1073 mVA' 

zdisp 

1 .96E-03 inches 

1.99E-03 inches 

ydisp 

3.33E-03 inches 

3.86E-03 inches 

max_stress 

73,457 psi 

74,720 psi 


Concluding Remarks 


A proposed method tor pertorming a structural 
optimization of a force balance using a computational 
experiment design has been presented. An example 
optimization ol a lorce balance measurement section 
has been illustrated demonstrating the feasibility of 
using this new method. 

The powerful insights that are provided to the 
balance design engineer using graphical optimization 
analysis cannot be over emphasized. This graphical 
optimization tool enhances the intuition and experience 
of the designer regarding design tradeofls and limiting 
boundaries created by the response criteria. Also, this 
type of analysis provides the balance designer with vital 
information regarding the complex, and often times 


subtle, underlying relationships between the design 
variables and the responses metrics. 

There exists considerable interest in an automated 
balance design system. The proposed method is well 
suited for an automated system by relying on objective 
design decisions. In order to automate this method the 
experimental design and analysis software would need 
to be integrated with the solid modeling and finite 
element analysis software. While there are significant 
benefits to an automated system, the cost of software 
integration would need to be Justified. 

The proposed method can lead to the optimization 
of current geometrical structural designs and can also 
be used to generate advanced geometry. The technique 
can be extended to include other structural design 
requirements such as dynamic response, which would 
combine the response metrics ot stillness and mass. 
Non-traditional response metrics could also be 
considered including manufacturing time and cost. 
While this paper focused on a specific structural 
optimization application, similar response metrics are 
commonly found in a broad array of structural 
optimization applications. 
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